Interaction of a CO molecule with a Pt 
monoatomic chain: the top geometry 
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Abstract. Recent experiments showed that the conductance of Pt nanocontacts and nanowires is 
measurably reduced by adsorption of CO. We present DFT calculations of the electronic structure 
and ballistic conductance of a Pt monoatomic chain and a CO molecule adsorbed in an on-top 
position. We find that the main electronic molecule-chain interaction occurs via the 5(7 and 2n* 
orbitals of the molecule, involved in a donation/back-donation process similar to that of CO on 
transition-metal surfaces. The ideal ballistic conductance of the monoatomic chain undergoes a 
moderate reduction by about 1 .0 Go (from 4 Go to 3 . 1 Go) upon adsorption of CO. By repeating all 
calculations with and without spin-orbit coupling, no substantial spin-orbit induced change emerges 
either in the chain-molecule interaction mechanism or in the conductance. 
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g INTRODUCTION 



The conductance of metallic nanocontacts produced by scanning tunneling microscopy 
(STM) or by mechanically controllable break-junctions (MCBJ) is routinely measured 
and reported in conductance histograms [1]. The nanocontacts are much shorter in size 
than the electron mean free path so that the conductance is ballistic and quantized, rather 
than diffusive and non quantized. Quantization is most readily observed in the form 
of peaks in conductance histograms obtained by superposing the conductance traces 
of many break junctions events. In some metals, such as gold and platinum, short tip- 
suspended monoatomic chains can be formed at the nanocontact [2-4]; in these cases 
the main histogram peak with the lowest conductance is due to the atomic chain. This 
peak is usually found near 1 Go (Go = 2e 2 /h) for clean Au nanowires [1,2] and at about 
1.5 — 2.0 Go in Pt, where 5d electrons contribute to the conductance [5, 6]. Accord- 
ing to Landauer's theory, ballistic conductance is determined by electron transmission 
and reflection at the contact. The presence of small adsorbates on the nanowire or in the 
nanocontact may naturally change electron transmission and reflection and alter substan- 
tially the conductance. In the case of Pt nanocontacts this effect has been demonstrated 
experimentally for H2 [7] and for CO [5, 6]. For instance, in presence of CO the con- 
ductance histogram of Pt is modified to exhibit two peaks at about 0.5 Go and 1.0 Go 
[5, 6]. Only the lower peak near 0.5 Go has been discussed theoretically [8], while that at 
1 Go has not. At the chemical level, in the experimental work of Kiguchi et al. [6] it was 
suggested that the relative strength of CO adsorption on Au, Cu, Pt and Ni nanocontacts 
could be explained by the Bly holder model [9], but the binding between a CO molecule 



and a Pt nanowire has not been addressed in this respect. Our calculations presently fill 
this gap. Finally, the effect of spin-orbit (SO) coupling on the molecule-chain binding 
and on the nanocontact ballistic conductance has not been taken into account in previous 
calculations, even though it is known that SO is so important to change significantly the 
electronic band structure of Pt monoatomic wires [10, 11]. 

In this work we study the interaction between CO and a Pt nanowire using a very 
simplified model: we consider a straight, isolated, Pt atomic chain with one molecule of 
CO adsorbed on-top (see the inset in Fig. 2). Calculations use a periodically repeated 
supercell containing iVpt atoms of Pt and the molecule with its axis perpendicular to 
the chain 1 . We perform electronic structure calculations in the framework of density 
functional theory (DFT) [12] within the local density approximation (LDA) for the ex- 
change and correlation energy using the PZ functional [13] (the optimized distances 
and the chemisorption energies have been calculated also with the gradient corrected 
PBE functional [14]). The role of SO coupling is investigated by comparing the results 
obtained with scalar-relativistic (SR) and fully-relativistic (FR) ultrasoft pseudo poten- 
tials 2 (US-PP) [16], the latter including the SO effect [15]. The ballistic conductance 
has been calculated with the Landauer-Biittiker formula, Go = 2e 2 /h T(Ep), evaluating 
the transmission T (obtained with the method described in Refs. 11, 17, 18) at the Fermi 
energy £p- All calculations have been performed with the computer codes contained in 
the quantum-espresso package [19]. 

ELECTRONIC STRUCTURE AND BALLISTIC CONDUCTANCE 

We optimized the platinum-carbon distance dp t -c an d the intermolecular distance dc-o, 
while keeping the Pt atoms aligned at their theoretical equilibrium distance in the iso- 
lated chain (2.34 A with LDA, 2.39 A with PBE). Within LDA the optimal geometrical 
parameters are d Pt - C = 1.82 A (1.81 A) and d C -o = 1.14 A (1.14 A) in the SR (FR) 
case, while the SR-PBE calculation gives slightly longer distances: dp t -c = 1.84 A and 
dc-o = 1.15 A. The adsorbed C-0 distance is slightly longer than the equilibrium value 
calculated for the isolated molecule, which is 1.13 A with LDA and 1.14 A with PBE. 
The chemisorption energy E c ^ em is calculated as the difference between the energy of the 
optimized geometry and the sum of the energies of the isolated wire and of the isolated 
molecule. With PBE, which usually gives realistic estimates, we get E chem = —1.4 eV, 
while LDA, which is known to overbind, gives larger values, E c ^ em = — 1 .9 eV in the SR 
case and E c ^ em = —2.0 eV in the FR case. Our results confirm that the interaction be- 
tween the CO molecule and the Pt chain can be rationalized by the Blyholder model (and 
its improvements), used to explain the chemical bond between CO and transition metal 
surfaces [9, 20]. The CO-metal bond is characterized by a balance between donation 



We verified the convergence of our results with respect to Np t and to the other parameters of the supercell. 
In particular the PDOS and the transmission have been calculated with iVp t = 17. 

2 Within LDA, the Pt pseudopotential is that given in Ref. 15, while for C and O we generated new US-PP. 
With PBE we generated new SR US-PP for all the atoms (the parameters will be reported elsewhere). We 
adopt a plane wave basis set with the following cut-off for the kinetic energy of the wavefunctions (charge 
density): SR -LDA, 29 (300) Ry; FR-LDA, 32 (300) Ry; SR-PBE 32 (320) Ry. 




FIGURE 1. Left: sum of the PDOS on the atomic orbitals of C (solid lines) and O (dashed lines) in the 
SR and FR cases. Right: SR-KS eigenstates associated to some of the peaks in the SR-PDOS. 



of 5(7 electrons from the molecule to the metal and back-donation from the Pt rf-bands 
to the empty 2n* molecular orbital, which gets partially filled by electrons. This can 
be argued by looking at the density of states projected on valence orbitals of C and O 
(PDOS), shown in Fig. 1 for the LDA case. The nature of the peaks in the SR PDOS has 
been identified through symmetry arguments or by directly visualizing the Kohn-Sham 
(KS) orbitals. In addition to the sharp PDOS peaks due to CO valence orbitals (3a, not 
shown since much lower in energy, 4a, 5a, Ik below the Fermi energy and In* 
above E?), we find peaks which are not present in the PDOS of the isolated CO and thus 
are due to the interaction with Pt. In particular, the peak at 1.2 eV corresponds to a new 
state 5o a (shown in Fig. 1) due to antibonding hybridization between the 5a orbital and 
the s and d bands of Pt; it lies above E? and thus gives rise to electron donation. The two 
sharp peaks at —4.9 eV and —3.5 eV, instead, are due to bonding hybridization between 
the 2n* orbitals and the d bands of Pt: since these new states lie below £p they are filled 
with metal electrons by back-donation. We note here that each of the % orbitals is split 
into an even state (27z£ e and 2k* & , both reported in Fig. 1, and l7r e ) and an odd state 
(l7T , 2n£ and 2n* Q ) with respect to the mirror plane which contains the chain and the 
molecule. We finally point out that, as in the case of CO adsorbed on surfaces, the 5a& 
level has moved below the In level, and that there is the same mixing between 4a and 
5a and between lit and 2n* (see Ref. 20). Comparing the position of the peaks in the 
SR PDOS with those in the FR PDOS, we can argue that this model of the interaction is 
valid also in presence of SO coupling. In fact, although the exact position of the peaks 
and their intensity may vary, we can recognize in the FR PDOS the same pattern of 
levels identified in the SR PDOS. This is due to the fact that the strength of the CO-Pt 
interaction is much larger than the effects of SO coupling on this binding. 

Fig. 2 shows the calculated transmission curve of the ideal (junction-free) Pt chain 
with the adsorbed CO molecule, calculated in the SR and FR cases within LDA. The 
calculated SR (FR) conductance is about 3.1 Go (3.0 Go), reduced by about 1.0 Go 
relative to the isolated Pt monoatomic chain, whose theoretical conductance is G = 4 Go 
(see Ref. 11). Below E?, between —1.6 eV and —0.4 eV we find a larger reduction 
in transmission with respect to the isolated chain, because the Pt states having orbital 
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FIGURE 2. Electron transmission as a function of the energy. The SR (FR) trasmission Tsr (7fr) is 
displayed with solid lines (open circles), while the transmission of an ideal monoatomic chain Nsr (Nfr), 
is shown with dashed (short dashed) lines. Inset: geometry and optimized distances. 



angular momentum along the chain m = ±2 are heavily perturbed in the vicinity of 
CO. Above Ep, we find a notch in the transmission curve at about 1.2 eV, precisely 
the energy where the new 5o a state is formed. The FR transmission agrees quite well 
with the SR one near E?, since the band with total angular momentum along the chain 
mj = ±5/2, which is raised in energy by the SO effect, is almost completely blocked 
by the interaction. In fact this band derives from the m = ±2 band, which is not well 
transmitted in the SR case. The largest differences between SR and FR transmission lie 
in the energy ranges —3.4 eV ^ E ^ —2.6 eV and — 1.7 eV ^ E ^ — 1.4 eV (where the 
different hybridization of the FR bands depletes the number of available channels), but 
they are far away from the Fermi level and hence they do not produce, in this system, a 
measurable effect on the conductance. 

REFERENCES 

1. N. Agra'it, A. L. Yeyati, and J. M. van Ruitenbeek, Phys. Rep. 377, 81 (2003). 

2. H. Ohnishi et al., Nature (London) 395, 780 (1998). A. I. Yanson et al., ibid. 395, 783 (1998). 

3. R. H. M. Smit et al., Phys. Rev. Lett. 87, 266102 (2001). 

4. V. Rodrigues et al., Phys. Rev. Lett. 91, 096801 (2003). 

5. C. Untiedt et al., Phys. Rev. B 69, 081401(R) (2004). 

6. M. Kiguchi, D. Djukic, and J. M. van Ruitenbeek, Nanotechnology 18, 035205 (2007). 

7. R. H. M. Smit et al., Nature (London) 419, 906-909 (2002). 

8. M. Strange, K. S. Thygesen, and K. W. Jacobsen, Phys. Rev. B 73, 125424 (2006). 

9. G. Blyholder, /. Phys. Chem. 68, 2772 (1964). 

10. A. Delin, and E. Tosatti, Phys. Rev. B 68, 144434 (2003). 

11. A. Dal Corso, A. Smogunov, and E. Tosatti, Phys. Rev. B 74, 045429 (2006). 

12. P. Hohenberg, and W. Kohn, Phys. Rev. 136, 864 (1964). 

13. J. Perdew, and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

14. J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865-3868 (1996). 

15. A. Dal Corso, and A. Mosca Conte, Phys. Rev. B 71, 1 15106 (2005). 

16. D. Vanderbilt, Phys. Rev. B 41, 7892(R) (1990). 

17. H. J. Choi, and J. Ihm, Phys. Rev. B 59, 2267 (1998). 

18. A. Smogunov, A. Dal Corso, and E. Tosatti, Phys. Rev. B 70, 045417 (2004). 

19. S. Baroni et al., website: http : / / www . quantum-espresso . org. 

20. A. Fohlisch et al., J. Chem. Phys. 112, 1946 (2000). 



